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Abstract 

We report the results of a detailed calculation of nucleon weak interactions rel¬ 
evant for the neutron to proton density ratio at the onset of primordial nucleosyn¬ 
thesis. Radiative electromagnetic corrections, finite nucleon mass terms, thermal 
radiative effects on weak processes and on neutrino temperature are taken into ac¬ 
count to reduce the theoretical uncertainty on n p rates to 1%. This translates 
into a sensitivity in mass fraction Yp prediction up to 10“'^. We find a positive 
total correction to the Born prediction 6Yp ~ 0.004. 
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1 Introduction 


Primordial Big Bang Nucleosynthesis (BBN) [Q, |^, represents one of the greatest success 
of the hot Big Bang Model and is among the main achievements of the interplay of particle 
physics and cosmology. The predicted abundances of all light nuclei produced in the early 
stage of the life of the universe, only depend on the baryon to photon density ratio rj, and 
a single value for rj ~ 10“^° may accomodate all observed order of magnitude for D, 
^He and '^Li data. In the last few years, due to the increasing precision in observations 0], 
an even more exciting possibility can be foreseen, namely that BBN may tell us new pieces 
of information on both fundamental physics beyond the Standard Model and our present 
understanding of the evolution of the universe. In this respect it is already remarkable 
that BBN has been able to strongly constraint the number of effective neutrino degrees of 
freedom. This perspective, however, the fact that BBN is now entering in his maturity and 
precision era, demands similar improvements in the precision of the theoretical analysis, 
in order to reduce as much as possible all uncertainties on the predictions. An increasing 
precision in the measurements of ‘^He mass fraction at the level of 10“"^ 0 requires, for 
example, a reliability of proton-neutron conversion rates at the level of one percent at 
least, and a control at the same level of precision of all other effects which are relevant 
for the neutron to proton ratio at the freeze-out. 

BBN took place at temperatures in the range O.l-i-10 MeV and can be seen as the hnal 
outcome of two stages, the decoupling of the weak interactions which keep neutrons and 
protons in chemical equilibrium and, shortly after, the onset of nuclear reactions which 
start producing light nuclear species. At temperature T >> 1 MeV the neutron to proton 
ratio is order unity and starts decreasing exponentially when the temperature reaches the 
value of their mass difference. As the temperature decreases, however, weak interactions 
are no longer fast enough to maintain equilibrium and a substantial hnal neutron fraction 
survives down to the phase of nucleosynthesis. All neutrons become practically bound 
in ^He nuclei, due to the high binding energy per nucleon, which therefore represent the 
dominant products of BBN, with an abundance in mass Yp close to 25 %. While Yp is 
weakly depending on the baryon to photon ratio rj, it is instead strongly inhuenced by 
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the ratio n/p aX the freeze-out, so it is particularly relevant to have a good theoretical 
estimate of this quantity. 

The nuclear reaction chain starts with deuterium formation Yiap+n —> D+'^ when the 
inverse rate for Deuterium dissociation starts falling. This happens when the temperature 
has decreased below the tiny D binding energy Eb = 2.2 MeV and, due to the large 
photon to baryon density ratio, when even the tail of photon Bose distribution becomes 
ineffective in D photodissociation, i.e. 77 “^ exp {—Eb/T) ~ 1, which gives T ~ 0.1 MeV. 
After Deuterium, the nuclear chain starts, leading to sensible production of nuclei up to 
Li (for a recent review see p). 

The more recent data on ^He mass fraction and D abundance are still controversial 
since there are two different sets of results mutually incompatible 

Yp = 0.234±0.0054 , D/H = (1.9±0.4)-10-^ ; 

Yp = 0.243±0.003 , D/H = (3.40±0.25)-10“® ; (1.1) 


which show the well known fact that large Helium mass fractions requires low values for 
Deuterium abundance and viceversa. It is possible that forthcoming measurements from 
high-redshift, low-metallicity QSO, or better understanding of the present data will clarify 
the situation , either in solving the controversy or pointing out a difficulty of standard 
BBN scenario. The fact which however is emerging from the above results is that Helium 
data are now reaching a precision of one per mille, requiring a similar effort in reducing 
the theoretical uncertainties on elementary n ^ p processes, as well as in evaluating the 
effect of uncertainties in the experimental nuclear reaction rates |^. 

The now classical numerical codes which are used to deduce nuclei abundances 0, 
have been already modihed by several authors, in order to include corrections to the tree 
level n ^ p transition rates. These corrections have distinct nature and, listed as they 
were historically introduced, are the following. 

First of all, when accuracy of the order of one percent is required, all Born amplitudes 
for n <->■ p processes should be corrected for order a radiative corrections. These effects 
have been extensively studied in literature | 


111 and can be classihed in outer factors. 


involving the nucleon as a whole, and inner ones, which instead depend on the details 
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of nucleon internal structure. Their inclusion in the /3-decay rate reduces the tree level 
estimate of neutron lifetime ~ 961 sec to the value ~ 894 sec, which is closer to the 
experimental one = 886.7±1.9 sec H- A similar correction, of the order of 10 %, is 
expected to affect the rates of all elementary processes relevant for BBN 
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Other small effects are actually expected at higher order, since the estimated value of is 
compatible with the experimental value at A-a level only. These additional contributions 
are usually taken into account by eliminating the coupling in front of the reaction rates 
(a)-(/) in favour of 

It was quite soon realized [|^ |l^ that since all reactions ( p..2|) take place in a thermal 
bath of electron, positron, neutrinos, antineutrinos and photons, hnite density and tem¬ 
perature radiative effects should be included as well. They may affect the Born rates for 
factors of the order aT/rrie, with rUe the electron mass, which in the temperature interval 
relevant for BBN may be as large as few percent. Calculations of these effects, using 
the real time formalism of hnite temperature held theory [^, has been done by several 
authors [T^, [0-[13]. Though they agree on the magnitude of the induced ehect on 
"^He abundance, there are however striking diherences in the sign of the corrections. 

Thermal ehects play also a role in modifying neutrino temperature after decoupling. 
The diherence between neutrino and photon temperatures originates from the entropy 
release to photon in e^-e“ pair annihilation which, roughly speaking, takes place at tem¬ 
perature of the order of the electron mass rUg, shortly after neutrino decoupling. The 
interaction with the surrounding medium gives an extra plasma mass term to e’''-e“, 
of the order of aT^/rrie, which slightly changes the neutrino to photon temperature ra¬ 
tio. Actually it has been shown that there is a residual small coupling of neutrinos at 
the electron-positron annihilation phase, which introduces non equilibrium effects in the 
neutrino energy distribution This last effect introduces corrections to ^He abun¬ 

dance of the order of (1 -i- 2)-10“^. 
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Finally all Born amplitudes for the processes in ( |1.2|) should be also corrected for 
nucleon hnite mass effects. This has been hrst considered in . They affect the allowed 
phase space as well as the weak amplitudes, which should now include the contribution 
of nucleon weak magnetism. Initial nucleons with hnite mass will also have a thermal 
distribution in the comoving frame, over which it is necessary to average the transition 
probabilities. 

In this paper we have considered in detail all the above effects, with the only exception 
of non equilibrium distortion of neutrino distribution, whose contribution to ^He mass 
fraction is however very small. In particular: 

i) all radiative corrections have been included at order a; 


ii) thermal radiative effects have been recalculated; 


hi) all hnite nucleon mass corrections at order l/Af/v are considered. 

The motivation for such an analysis which, as we have stressed, is certainly not new, is 
twofold. On one side for some of these ehects diherent results have been obtained by 
several authors. Furthermore the precision expected for ^He mass fraction measurements 
requires a careful scrutiny of the net ehect of all these one percent level corrections. 

We have organized the paper as follows. In section 2 we report the customary ex¬ 


pressions for Born rates for the processes (|1.2|) . Section 3 is devoted to radiative zero 
temperature corrections. In section 4 we discuss hnite nucleon mass corrections, while in 
section 5 we report all our results for thermal radiative ehects, evaluated in the real time 
formalism. In section 6 we present our results and discuss the expected correction to ^He 
mass fraction. Finally in section 7 we report our conclusions and outlook. 


2 The Born rates for n ^ p reactions 

At the epoch of BBN, 0.1 MeV < T < 10 MeV, the universe was hlled by a primordial 
plasma of nucleons, electrons, positrons, neutrinos, antineutrinos and photons, which were 
initially in thermal equilibrium. The initial conditions for the primordial nucleosynthe¬ 
sis are settled out by the ratio of the relative abundances of neutrons and protons. For 
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temperature larger than 1 MeV this ratio is strongly hxed by thermodynamical equilib¬ 
rium, but as T continues decreasing the necessary conditions to sustain the equilibrium 
are no longer satished. In this regime, to obtain the neutron-proton abundances, it is 
necessary to solve the Boltzmann equation, where the rates for the charged-current weak 


interactions transforming n p (|1.2|) are involved. Shortly after, when the temperature 
drops down to 0.1-^ 0.3 MeV, the nuclear reactions become strongly active and transform 
neutrons mainly in and only in a small fraction in the other light elements. 

To obtain the rates per nucleon for the n p processes, let us consider a generic 
transition i —>■ f in (|1.2|) . The differential reaction rate can be written in the form 


dT{t^f) = V |M(z ^ /)P TT . F, n , (1 - Ff) 

^ ^ V(2vr)32po. (27r)32po/^ 


spins 


X 




( 2 . 1 ) 


where \M{i —/)p is the squared matrix element, to be summed over all spin degrees of 
freedom, the four-momentum of j-th particle is pj = {poj, pj), and Fj denotes the phase 
space density of the j-th particle. 

Since the temperatures involved during the BBN are less than 10 MeV, in hrst ap¬ 
proximation one can neglect in ( p.l|) the kinetic energy of nucleons with respect to their 
mass. Moreover the Pauli blocking factor for the hnal nucleon is completely negligible, 
since nucleons are strongly nonrelativistic and < exp(—M atj^/T) < exp(—100), where 
Nf denotes the hnal nucleon. In this limit, the differential reaction rate per incident 
nucleon {Ni) becomes 


duj 




— Ott \M{i —» /)p y-j- 


d?pi 


spins 


F^ n 




(27r)32poi 


(1 - Ff) 


X 6 I MjVi — MjVj. + '^ Poi — X! Pof 

i^Ni f^Nj 


( 2 . 2 ) 


Hereafter, we will denote with uj{i —>■ /) the reaction rate per incident nucleon for the 
process i ^ f and the index B reminds that the rates are computed at tree level and in 
the inhnite mass limit for nucleons. We will refer to this as the Born approximation. 

Let us start considering as an example the reaction (a) in (|1.2|) . At tree level, the 
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corresponding V — A amplitude for the diagram of Figure 1 is 


G 

M{ue +n^ e~ +p) = - C^75)Mn(g') Ue{p')l^{l - l5)uu{.q), (2.3) 

where Gf is the Fermi coupling constant and Cv, Ca are the vector and axial coupling of 


the nucleon. In Born approximation we have for (p.3|) 


E 


|M|2 


. 8MnM„ 

spins P 




(cl + 3G^) Po^o + (cl - Cl^ p' 


■q 


(2.4) 


Thus, by substituting (|2.4|) in (p.2|) and after integrating we get 

Gj,{Cl + 3Cl) r 


+ n^ e +p) = 


271 ^ 


d\p\ IpT % 0(9o) Fy{qo) [1 - Fe(po)] , (2.5) 


where the integration limits are imposed by the ©-function, go > 0. In the present case we 
have go = —{Fin — FIp) + p^ = —A + Pq. The Fermi statistical distributions for electrons 
and neutrinos in the comoving frame are [] 


Fe(p') = + 1 


n -1 


Fu{qo) = + 1 


n -1 


( 2 . 6 ) 


with (3 = 1/T and jSy = l/T^. 

Due to crossing symmetry, for the other hve processes in ( p..2|) the form of I^spins |Tfp 
in (p.2|) remains unchanged, so the hnal result can be obtained from (|2.5|) by simply 
replacing go, determined by the energy conservation for each reaction, and the thermal 
factors, which depend on the initial and final lepton states for each process. In Table 1, 
using the same notation adopted for z/g -|- n —»• e“ -|- p, we report the values for go and the 


thermal factors for each reaction in (|1.2|) , while in Figure 2 we show the Born rates for 
processes ( |1.2| ) versus photon temperature. 

The electron neutrino temperature behaviour after decoupling is shortly reviewed in Ap¬ 
pendix A. The result quoted there for the ratio Ty/T will be used all through our analysis. 
Non-equilibrium effects on neutrino distributions due to residual coupling during the e^- 
e~ annihilation phase, and their implications on ^He abundance will be shortly discussed 
in section 6. 


^ We are assuming vanishing chemical potential for neutrinos. See for a review on the effects induced 
by neutrino asymmetry and constraints on chemical potentials. 
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Table 1. 
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3 Radiative electromagnetic corrections: the neutron 
lifetime 


In order to measure the level of validity of Born approximation in the evaluation of rates 
(PD , one can use the corresponding prediction for neutron lifetime. By using (|2.5| ) and 
Table 1, in the limit of vanishing density one gets the simple expression 

-ml I dee (e — 


, ^ Gl {Cl + 

27r3 ^ 


- -)' - 1 ) 


(3,1) 


Thus, inserting the value of Gp, very well known from the measurement of muon decay 


T^, Cy = Vud = 0.9751±0.0006 and CaICv = 1.2601T0.0025, deduced from the study 


of the decay product angular distribution in neutron decay |^, we obtain the value 
Tn — 961 sec which has to be compared with the experimental result = 886.7 ± 1.9 sec 
12]. As a consequence, a correction factor to (H) of about 8% is needed. Note that 


similar corrections are expected to affect the other reactions in (|1.2|) as well, since they 
are mediated by the same interaction hamiltonian. 

Let us consider order a QED corrections to neutron lifetime {a is the fine structure 
constant); these can be separated into otifer corrections, involving the nucleon as a whole, 
and inner corrections, depending on nucleon structure. Obviously, the inner corrections 
are sensible to the details of the strong interactions inside the nucleon, while it can be 
shown P] that outer corrections, at least up to terms of order a, take a general form. 



























In terms of Feynman diagrams the corrections at order a to nentron lifetime are 
obtained from the pure V — A interaction ( |2.3|) with the exchange of a photon between the 
four fermions involved. Note that nucleons interact electromagnetically through both their 
charge and magnetic moment, even though the interactions involving magnetic moments 
are suppressed by the inverse of the nucleon mass. 

The outer correction to the differential decay rate can be written as the multiplicative 
factor 

(3.2) 

where g{pQ, go) is a function of electron and neutrino energy, whose explicit expression 
can be found in eq. (206) of Ref. [Q. This function describes, at order a, the deviations 
in the allowed electron spectrum arising from the radiative corrections^ It increases the 
decay probability for neutron /d-decay (and therefore decreases neutron lifetime) by about 
1.6 %. 


a 


^ + ^9{Po,go) 


Inner corrections are much more difficult to handle since they strictly depend on 
nucleon structure. In general, one can follow two different approaches. One possibility is 
to directly consider radiative corrections to the effective nucleon weak current Mp 7 ^(l — 
{Ca!C v)'^^)Un Alternatively, one can study corrections to quark currents ? 7 /i(l— 75 ) 5 ' 

and then translate the quark-based description into the hadronic one 0. Here, we adopt 
the second point of view, and report the results obtained by Marciano and Sirlin P 


o. (. , Mz , Mp 
— \ + 2C + A. 


27r 




Ma 


(3.3) 


j.p ) 

The hrst term represents the dominant model-independent short-distance contribution. 
The second and third terms are axial-current induced contributions, where Ma is a low en¬ 
ergy cutoff applied to the short-distance part of the 7 IF box diagram and 2C is the remain¬ 
ing long-distance (low energy) correction. These terms depend on the details of the strong 
interaction structure and are the main sources of uncertainty on the radiative corrections. 
The allowed range for the cutoff Ma is 400-71600 MeV, while C = 0.798 Ca {l+Pp + Pn), 


At is customary to substitute the function g{pQ, go) in (^) with its mean value g obtained averaging 
over the allowed electr on sp ectrum. However, since we will consider not only neutron decay, but also the 
other five reactions in (H2), we will not use this approximation. 
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/^p + /^n — —0.12 being the nucleon isoscalar anomalous magnetic moment. Then we have 


a 


— In 
27r 


(in^ + 2C^ ~ 0.0012 ±0.0018 . (3.4) 

The last term in (O) is a perturbative QCD correction whose calculation is rather reliable 
and gives Ag ~ — 0.34. 

The largest correction comes from the hrst term in ( |3.3| ) and it seems appropriate to 
approximate the effects of higher orders by summing all leading-logarithmic corrections 


of the type a”ln"'(M^) via a renormalization group analysis. This has been done in 
the corrected differential rate then acquires the multiplicative factor 


GiPo^ Qo) = 


1 ± 


a 




Ma 


a{Mp 


In ^ ± 2^ ) H-[5((po, qo) ± Ag] 


27r 


S{Mp,Mz) 


(3.5) 


where a(/i) is the QED running coupling constant dehned in the MS scheme. The short- 
distance enhancement factor S{Mp, Mz) is given by |]l0 
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S{Mp,Mz) = 


a [me) 

a{Mp)^ 


aimr 


aime) 


_ 9 _ 
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a{mb) \^ 1 ® / a{Mw) \ / a{Mz) \ 

\ a{mt) ) \a{Mw)) 


17 


aim. 


. (3.6) 


In the MS scheme the running coupling constant at different scales (a ^(M^), a ^{M]y), 
a-\mb), a-\mr), a“^(mc), a-^{Mp)) is (127.90, 127.94, 132.08, 133.27, 133.62, 133.90) 
leading to ^(Mp, M^) ~ 1.0225. 

The radiative corrections just considered, which are described by the factor in (|3.5|) , 
are not, however, the only ones to be taken into account for the neutron decay. They 
contribute for about 4 % to the lifetime, against the required global correction of about 
8%. The remaining leading corrections are usually viewed as corrections to the phase 
space factor |^, namely to the differential rate in (|3.1| ), since they are electron and 
neutrino energy dependent. 

Let us hrst examine the correction due to the distortion of the outgoing electron wave 
function by the Coulomb held of the proton. This can be calculated by solving the Dirac 
equation for an electron under the inhuence of a spatially hnite proton charge distribution 
of radius i? ~ 1 /m; the wave function is then evaluated at the centre of the proton. The 
correction factor is then given by [pT 


^{Po)^iPo) - (1 


P' 


m: 


1 — aRp'n 1 -I-9 

V Vo^ 


(3.7) 


^The expression reported in H has been updated to account for the measured top quark mass. 
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The first term is usually denoted as the Fermi function for the Coulomb scattering. 

Another small correction comes from the observation that neither the electron wave 
function, evaluated at the centre of the proton through J^{Pq) jC,{pq), nor the anti-neutrino 
one are constant through the proton volume. Thus, the decay rate has to be computed 
by an appropriate convolution of the electron, anti-neutrino and proton wave functions 


through the proton volume. This leads to the following correction factor |11 


C{Po,qo) ~ 1 + Co -h ^ + C 2 P 0 + C 3 P 0 

Po 


/2 
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C 2 = 
C3 = 
B = 


mill 
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(Po + Qo? 
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aBR{p'Q + go) 


2B 2/ / N 
— R m,(po + go) 




, 3 (Po + 9o) - “ 


-(3 - B)R^ 


Cl 


Cl 


(3.8) 

(3.9) 

(3.10) 

(3.11) 

(3.12) 

(3.13) 


(C'i + 3 C 1 ) 

Finally, there are also tiny corrections related to the fact that the outgoing proton is not 
at rest. An exhaustive discussion of hnite nucleon mass corrections will be given in the 
following. Here, we only point out that since proton recoils, its Coulomb field is the one 


of a moving source. The correction associated to this effect is contained in the factor 

,2 / 


Q(Po,9o) ^ 1 


airniz 


1 + B 


Qo 


Mp\p'\ V ’ ^PoJ ' 

We eventually get for the neutron /3-decay differential rate, corrected at order a 


(3.14) 


duj(n -^e + +p) = e(p!i,«o).?'(p!i)2:(p!,)C(pi,(/o) Q(p'o,qo) di^B(n ^ e +v^+p), 


(3.15) 


and thus the neutron lifetime reads 

Gl {Cl + 3 ^ 1 ) 


= 


27r^ 


dp'oP'oW\ {p'o - A)^ 


X G{p'q, Qo) d^iPo) ^(Po) qo) Q{Po, go) 


(3.16) 


Concerning the relevance of the different corrections just considered, a comment is 
in turn. For seek of completeness we have reported an exhaustive description of all 
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corrections at order a, but as clearly appears from the explicit expressions of (|3.5| )-( |3.14|) , 
the main contributions come from ( ^.51) and J^(Po) term in (|3.7|) . All other contributions, 
contained in C{pq), C{pq, go) and Q{Pq, qo), are in fact much smaller and thus they can be 
safely neglected, since they are suppressed by a factors of the order A/Aqcd- Actually 
the contributions coming from radiative and Coulomb effects are even larger, or of the 
same order of magnitude, than these terms, so they should be included for consistency if 
C{pq), C{pq, go) and Q{Pq, go) are taken into account [j^. For the level of accuracy of our 
analysis it will be sufficient to include J^{pq) and G{pq, go) at order a only. 

Evaluating numerically the integral over the phase space we then obtain for the neutron 
lifetime the value Q = 893.9 sec, which is now compatible in about A-a with the 
experimental value. The theoretical prediction for neutron lifetime can be further refined 
by considering, for example, higher order terms in a, like subleading corrections, or other 
small effects, such as residual average proton polarization due to parity non-conservation 
and so on. Some of these effects are briefly discussed in [|^], but for our purposes we do 
not consider them, since their small contributions would be mainly energy independent 
and then account only for an overall factor. In fact, assuming that such residual improving 
would eventually yield a theoretical prediction for neutron lifetime fully compatible with 
the experimental value, it is meaningful, as far as the aim of the present investigation 
is concerned, to rescale the differential rate of neutron /3-decay in (|3.16|) for the residual 
constant small factor This is a standard procedure, which allows to overcome 

the problem of a precise determination of the coupling constant for weak interactions 
involved in BBN, expressing this overall factor in terms of the experimental value of 
neutron lifetime. We stress, however, that it is worth rescaling the rate only after all 
known corrections have been included, since this increases the accuracy of the prediction. 
In this way the ansatz that the residual correction be an overall factor and not, as it is 
reasonable to expect, a function of leptons energies, may introduce errors less than 1%. 

The analysis of radiative corrections for /3-decay shows that their contribution to 


the rates for the processes in (|1.2|) , relevant for BBN, is expected to be as large as 
few percent. To reach accuracy of the order of one percent in rate evaluations it is 


^For completeness we have also included corrections coming from the finite nucleon mass, which are 
considered in section 4. However, these corrections are very small. 
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therefore necessary to correct Born rates for both radiative and Coulomb effects. This 
last contribution, however, is only present when both electron and proton are present 
in the initial or hnal state, since it can be viewed as the electromagnetic rescattering of 
the two charged particles. This results can be also obtained diagrammatically, using the 
traditional separation of electromagnetic corrections into order a (radiative) and order 


Za (Coulomb) contributions, though this separation for the elementary reactions in (|1.2|) 
is somehow arbitrary, since Z = \. In this way it is straightforward to show that no 
Coulomb corrections .F(pq) are present for the channels e"*" + n pg (see and 
references therein). 

Summarizing, zero temperature radiative corrections modify the Born elementary rates 
in the following way: 


i) the radiative energy dependent factor Qip'^^qo) reported in (|3.5|) is included in the 
differential rates for all processes {a)-{f); 

ii) the Coulomb factor (^3) introduced for processes (a), (6), (e) and (/); 


iii) all rates are rescaled by the factor 


1 + 5, = rf 


6 ^ = 0.008 


(3.17) 


representing an energy independent constant correction which, as we have discussed, 
should be included to reproduce the experimental value for neutron lifetime at zero 
density. 

4 Finite nucleon mass corrections 


In this section we consider the corrections to the Born rates for the n ^ p reactions 


(|1.2|) induced by relaxing the assumption of inhnite mass for nucleons. This is a necessary 
improvement since the hnite nucleon mass corrections, proportional to rUe/MN or T/Mn ~ 
a, are of the same order of magnitude of the radiative ones. 

The main effect ^9 of considering hnite values for the nucleon masses are to allow for 
new effective extra couplings, such as weak magnetism, and to change both the evaluation 
of the spin summed squared matrix element and the allowed phase space. Eq. (p.4|) is in 
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fact only valid for nucleons at rest. Furthermore, the fact that all rates should be averaged 
over the Maxwell-Boltzmann distribution for the initial nucleon, induces corrections of 
order T/Mjy. This kinematical correction will be considered separately. 


4.1 Corrections to the transition amplitude and phase space 


For dehniteness let us consider the process shown in Figure 1. At hrst order in 1/Mn, 
the transition amplitude can be written as 

Gf_ 


M = -^Up{p)O^Un{q')Ue{p')Y{^-l5)Uu{q) , 

where the effective nucleon-nucleon weak coupling is given in general by 


(4.1) 


= 'Jp{Cv - + i {p - q'Y 

Mn 


(4.2) 


The form factor /2 is the anomalous weak charged-current magnetic moment, and we have 
neglected both scalar coupling, related to CVC breaking, and pseudoscalar coupling of 
the nucleon, whose small value is hxed by low-energy PCAC theorems ||^, [^|. For the 
weak magnetism term / 2 , we have from CVC, /2 = Vud{Pp — Pn)/‘^ = 1-81 Vud- In general, 
the form-factors Ca and / 2 , are all dependent on the transferred four-momentum. 
However, at the energy scales relevant for the considered processes, this dependence can 
be safely neglected. 

In appendix B we have reported the spin-summed squared modulus of the transition 
amplitude (|4.1|) . The result is cast in a completely general form in order to allow for a 
straightforward generalization to all processes of ( |1.2|) . Specializing this result for the 
process Ue + n ^ e~ + p, and including pure radiative effects analyzed in the previous 
section we get the following expression for the rate 

1 + 


ujR{h'e + e +p) = 


mMtJo Slvp(pi.«o) 


spins 


X S(Po,qo) Fdqo) |1 - felPci)] 


(4.3) 


where 


Einf 


{Mn + qo){Ml - Ml + ml + 2qoMn) - 2qo^ 


2 Mn{Mn -|- 2go) 


(4.4) 
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(4.5) 


E — 

^sup 


i = 


{Mn + go)(^n ~ 


2 Mn{Mn + 2 go) 


M2 - M2 - m2 + 2goM„ 


1 2 


— AmiMf 


(4.6) 


Notice that the squared amplitude is evaluated in the neutron rest frame, and that the 
result is not averaged over the neutron velocity distribution. This effect will be considered 
in the following section. In Figure 3 we plot the zero-temperature radiative correction 
/S.UJR = ojr — ujb normalized to {a/Ti)ujB for the processes ( p..2D . Weak magnetism and 
hnite nucleon mass effect on phase space are included. The constant shift correction vrh ,-/a 
(see eq. (|3.17|) ) has been subtracted in order to show the pure radiative and Coulomb 
effects and hnite mass corrections. 

4.2 Kinematical corrections 

Apart from mass corrections to the squared matrix element and allowed phase space, there 
are indeed additional contributions to the reaction rates which are of a purely kinematical 
nature, namely the ones connected to the kinetic energy carried by the initial nucleon in 
the comoving frame. 

Let us consider the reaction + n ^ e~ + p and, with the same notation of the 
previous sections, its rate per initial nucleon 


UjiVe + n 


p) = — 




jFniq'o, 


d^p d^q d^p' 


(27r)32g' [J (27r)32po (27r)32go (27r)32p' 

X + q' - p - p') ^ |M|4 F^(go) [1 - F;(po)] , 


(4.7) 


spins 


where, as already mentioned, one can safely approximate 1 — Fpipo) ~ 1 and 

d3q/ 


Tlri, 2 


(2ir) 




(4.8) 


denotes the incident neutron number density. All quantities in ( [4.7| ) are evaluated in the 
comoving frame, where the plasma is at rest, which coincides with the initial nucleon 
rest frame only in the inhnite nucleon mass limit. Let us hrst note that the quantity in 
brace brackets is a Lorentz invariant so that it can be equally evaluated in the neutron 
rest frame. If we therefore denote with % and p'^ neutrino and electron energies in the 
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neutron rest frame and consider the change of variables in the integral ( |4.7| ) ( 5 'o,Po) —>■ 
{Qo,Po), the transformed expression in brace brackets is simply obtained by substituting 
<?o —^ Qo, Po Po- To evaluate how electron and neutrino thermal distribution factor 
Th((?o) [1 ~ ^e{Po)] transforms under this change of variable, a simple computation gives, 
up to 1/M^ terms 


Po - Po + 


Mr,. 


+ 


|q 


/|2 


2 M 2 


Po 


q'q 

* » + ^ + 


2 M 2 


% 


and thus one gets 

FMo) [1 - F.(pi)] = -R.(9 o) |1 - f.(pi)] + [F, (1 - F.)] 

where 

[1 - )] 


(4.9) 

(4.10) 

(4.11) 


FM) - [1 - FAPo)] 


|q 


2 M 2 T 


PoFeiPo) - 


|q 


2 M 2 T, 


M„,T 


go [1 - F^{qo)] - - 


Mr,. T,. 


' q^ p- - 
Mr.T . 


n IP 
2 


FeiP',) [1 - 2Fe(p(,)] 


q'-p'q'-q r /- m ,1 

“ F,,{qo)] + - 


q'q 


[1 - T;(?o)] [1 - 2F^(go)] 


(4.12) 


By substituting (|4.11D in ( |4.7|) it is easy to see that the hrst term on the r.h.s. of (|4.11|) 
reproduces the results of the previous section. Hence, the new contribution, proportional 
to 5 [Fy (1 — Fe)] of ( [4.12|) is an additional 1/M„ correction to the Born rate which van¬ 
ishes for |q'|/M„ —>■ 0 and is of a purely kinematical origin. It is due to the fact that 
the comoving and the neutron rest frames are not coinciding, unless neutron are taken 
inhnitely massive, and initial neutrons have a velocity distribution in the comoving frame, 
over which the reaction differential rates should be averaged. 


Using that in the nonrelativistic limit 


Fn(q'„) = 


exp 


'i) + il" 


~ exp 


T 


~ exp 


' T 


exp 


|q 


/|2 


2 M„T^ 


(4.13) 


the kinematical correction is 


+ p) = 


_ Mn 

e F 


d^q' 


|q' 


/|2 


Un J (27r)32g(, 


g 2 MnT 


d^p d^q d^p' 


(27r)32po (27r)32go (27r)32p' 


X {27r)U^^\q + q' -P-P') ^ |M|n d [F, (1 - F^)]. (4.14) 


spins 
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Since S[F^{1 — Fe)] is already a 1/M„ correction, one can compute the above integral 
assuming for the quantity in brace brackets the inhnite nucleon mass limit. According to 
this approximation, the dependence of the integrand function on q' comes through the 
expression (^4.12|) of 6 [F,^ (1 — F^)] only. 

Performing in (f4.14|) hrst the integral in and observing that 




exp 


|q 


2M„T 


and 


_ Mn 

e T 
n„. 


■ exp 


|q 


we get the result 

+ n ^ e~ +p) = 

x[l-Fe(Po)] 


d^q' 

(27r)32M„ 

GUCl + ^Cl) f T 


C 2 


2M„T 


q' = 0 , 


T 


(4.15) 


(4.16) 


27r^ 


2MnJ Jo 


d\p'\ IpI^^o ^iQo)FMo) 


-3 + ^ [1 - FJqo)] |1 - 2FJq„)] - ^F,(pi) [1 - 2F,{p'„)] 


Cl-Cl 


^ r) I /1 2 

2|p I 9o 


2^^(9o)l + 3j,Fe(Po)+ y c'j + scjy \3TT^pi 


2 FM) [1 - ^’.(9o)l 


. (4.17) 


The above expression for + n ^ e + p) cannot be extended to all reactions of 

( |1.2|) by simply using the substitution rules of Table 1. This is due to the different form 


of the statistical factors involved in the corresponding expressions of (|4.7|) . In Appendix 
C we have reported the corresponding correction terms for the other processes (&)-(/) in 
( |1.2|) . Using these expressions the contribution Auk can be numerically evaluated. We 
report the results in Figure 4. 


5 Thermal radiative corrections 


In the evaluation of the n ^ p rates, in addition to the pure electromagnetic radiative cor¬ 
rections, already analyzed in the previous sections, it is also necessary to include thermal 
radiative corrections 10.0- These are due to hnite temperature effects in the radiative 
processes of Figure 5, whose strength is expected to be of the order of a(T/me), and are 
related to the interactions of the in/out particles involved in the microscopic processes 
(a)-(/) with the surrounding plasma of pairs, r'e-W and photons. As in previous 

sections, for seek of brevity we will explicitly illustrate the computation of corrections 
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for the particular process z/g + u —>■ e~ + p. The generalization to all six processes is 
straightforward using Table 1. 

To compute such thermal effects it is useful to apply the standard real time formalism 
for hnite temperature quantum held theory |T^. In this approach propagators for 7 , 
e^, Ue and Vg appearing in Feynman graphs of Figure 5, take the following expression, 
containing a hnite density term as well 


i Seip') 
iSy{q) 


kP' — 1 

,/2 


9^'^ = - 




+ 27r 5{k^) B{ko) 




- me 


271 6{p' - ml) Fe(po) (/ + me) , 


27r 6{q‘^) Fy{qo) (f . 


(5.1) 

(5.2) 

(5.3) 


The additional contributions^ pick up real particle via the mass-shell h-functions and are 
proportional to the particle densities in the thermal bath. They take into account the 
role played by the medium in the single particle propagation. Note that a similar term 
for nucleons would be suppressed by a Boltzmann factor smaller than exp(—100) in the 
temperature range relevant for BBN, so it will be neglected in the following. Hereafter 
we will denote with thermal radiative corrections the hnite temperature ehects for all 
processes corresponding to the graphs in Figure 5, once the pure radiative contributions, 
considered in previous sections, have been subtracted. 

The diagrams of Figures 5(a), 5(b) provide corrections of order a{T/me) via the in¬ 
terference terms with the corresponding tree-level diagram. In particular, the diagram of 
Figure 5(a) is responsible for electron mass shift and wave function thermal renormaliza¬ 
tion terms. The diagram of Figure 5(b) provides the vertex correction via its interference 
with the tree-level process. Finally, in order to cure the infrared/collinear divergences 
present in the 5(a) and 5(b) contributions, it is also necessary to consider the squared 
amplitudes corresponding to the diagrams of Figures 5(c)-5(f), which describe photon 
emission and absorption processes. 


®We remind that we are assuming vanishing chemical potential for all species. 
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5.1 Electron mass shift and thermal wave function renormaliza¬ 
tion 

As already stated, the diagram of Figure 5(a) is the electron self-energy contribution to 
the Pe + n ^ e~ + p weak process and contributes to thermal radiative corrections via 
electron mass and wave function renormalization. The electron propagator reads 

1 


5.= 


— rrif. — S ’ 

where the self-energy S is given by 

r 

S(p') = idvra J ■ 

Hence, the thermal part of S is 


(6.4) 


(5.5) 




2771^ — '^' + /: ^ xn2\ 2me + 


ip' — kY — ml 


B{ko) 6{k^ 


{p' + ky 


F^E) 5{P - ml 


( 6 . 6 ) 


with E = Y |k| -I- ml- 

By choosing the direction of electron propagation along the 3-axis, the thermal self¬ 
energy takes the simple general form 


Tit — So -|- 7 ° Si -|- 7^ S2 , 


(5.7) 


where the explicit expressions of the Sj can be deduced by comparing (p.6|) with (|5.7|) . 

Differently from the usual self-energy contribution, the hnite temperature contribution 
contained in is not Lorentz invariant. This is due to the presence of the particle 
densities B{kQ) and Ff,{E), which take the explicit expressions reported in (|5.6|) only in 
the comoving frame. Actually, the lack of Lorentz covariance, due to the presence of 
the preferred frame of the medium, makes tricky the dehnition of a particle state, which 
is generally dehned as an irreducible representation of the Poincare group. All these 
difficulties originated a long debate in literature [jlS, 0, [|^-|^, concerning the correct 
way to evaluate the thermal wave-function renormalization, which is strongly affected 
by this ambiguity. Recently, a simple approach has been developed in which the 
dehnition of particle naturally comes out as the energy pole of the perturbed propagator. 
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This approach leads to the following expression for the renormalized projector on positive 
energy states 






1 


(% , Si\ 

r 1 

/ , P 

m —r 

P-7 + 

VIpI po) 

L me \ 


. (5.8) 


In the previous expression p'J^ = \/|p 1^ + These notations will 

be implicitly used hereafter. The hat on Sj and S' = dTii/dp^ reminds that they are 
evaluated on shell, i.e. for p'q = + rrig. Furthermore, the electron renormalized 


mass m!} results to be 


TTl^ — ?7le + Sq H-^Si + - - '-^2 

rrip nip 


P'o 


IP 


(5.9) 


A similar result is obtained for projector on negative energy states 0 . 

As it is clear from (|5.8|), at finite temperature the projector operator does not renor¬ 


malize multiplicatively. The presence of an additional term in the expression of A^ was 


already pointed out in literature |^, but the form of this term did not have the correct 
limit in the T = 0 covariant case. This is instead the case of (|5.8|) , as it can be simply 
verified. However, the additional term in (|5.8|) does not affect the electron wave-function 


renormalization contribution to rates for processes (| 1 . 2 |) in the nonrelativistic limit for 
nucleons, as can be seen by explicit calculation^. 

For brevity we only report the expressions])] for Sj involved in ( b.9|) 

MFe{E) 


50 

51 

5 2 


anip 

ttIp'I . 


d|k| 


6 |h 

iraT'^ 


^log 


E 

'po + Ip'P 


(log A — logH) 


,Po 


IP' 


+ 


a 
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3|P'| 
1 - 


2 |p'| ^ 


27 r|p'| Jo 

Vo + IpT ^ 
- Ip'L 


d|k| \k\Fe{E) {log A + logB) 


(5.10) 

(5.11) 
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2q; r= 
7 r|p'| Jo 


d\k\ 


\^fFe{E) 
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PqE + ml ^ . p'nE - ml ^ 

^ Mog A - ^ log B 


4|k||p'| 


4|k||p'| 


with 


A = 


p'qE + ml + \ 

|k| 

Ip'I 

^ _PoF-ml+ \ 

k| 

Ip'I 

PqE + ml — \ 

|k| 

Ip'I 

p'oE -ml- 

|k| 

Ip'I 


( 6 . 12 ) 


(5.13) 


®To evaluate thermal radiative corrections one can safely assume the approximation of infinite mass 
for nucleons. 

^We note that our result for E 2 differs for a numerical factor from the analogous one of Ref. |14|. 
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For the process e +p the Born graph yields the expression for the rate for nucleon 

given by ( 0 ). The additional contribution due to the electron mass renormalization Acum 
can be obtained by replacing rrie —>■ in ( ^TSf ) ; thus at order a we get 


AujM{i^e + n ^ e + p) = 


GUCl + 3Cl)a /•- 


X 


271^ TT 

^ , FejPo) _ 1 - F^jqo 

qo T T, 


/*oo 

/ d\p\\p'\ Qo Q{qo) Fu{qo) [l-Fe{po)] 

Jo 

2 Fe{E) 
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H—r 




Po F 


d|k| |k| 
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mt 


2 Ip' bo Jo 


dlkl lkl 


Fe{E) 

E 


(log A — log-B) 


(5.14) 


Note that the extension of the above result to all reactions of (1^) cannot be simply done 
using the substitution rule of Table 1, since the dependence on the renormalized mass 
comes through the statistical factors of the Born expression as well, which depend on the 


particular reaction. In Appendix C we have reported how (|5.14|) must be modihed to be 
applied to the other processes ( 6 )-(/). 

The contribution due to the electron wave function renormalization is obtained by 
using the projector (^.81) in the evaluation of the Born rate. This procedure, having 
subtracted the contribution due to the mass renormalization gives, at order a 

GUC^^ + 3Cl)a r 


+ +p) = 


2rp2 
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ttV 


7iv‘°S 


27r^ TT 

'p'o + Ip' 


. 6 |p'bb 
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2 |p'bb 2 o 


d|k| |k| 


-IP' 

Fe{E) 
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+ 2 


/*oo 

/ c^Ip'I Ip'I ql 0 (9o) F^i^qo) (1 - Fe(po)) 

Jo 

2 Fe{E) 1 


d|k| |k| 
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IP' 


[{E + Po) log A+{E- p'q) log B] 


(5.15) 


In the expression ( |5.15|) we have also subtracted an infrared divergent term which is 
cancelled by an opposite term provided by the photon absorption and emission rate. 
The collinear divergence, still present in (|5.15|) for |p'| = |k|, together with a similar 
contribution coming from the vertex correction is also compensated by bremsstrahlung 
diagrams. 

Differently from (^.14|) , the extension of (|5.15|) to other processes can be straightfor¬ 
wardly obtained by using Table 1. The results for Aum and AojWi normalized to {a/Ti)ojB 
versus photon temperature are reported in Figures 6 and 7, respectively. 
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5.2 Vertex Corrections 


The vertex correction is provided, at order a, by the interference term between the dia¬ 
gram of Figure 5(b) and the tree amplitude of Figure 1, leading toQ 

Auvii^e + n^ e~ + p) = ^ _ f d\p\\p'\f d|k| |k| 0 (go) 

tt Jo Jo 

FeiE) 


xFuiqo) [1-F;(Po)] 


E 


\ogA + — 


E 


log 5 — 


2|k| 

Ip'I 


Ip'I 

2 


k 

2 


(5.16) 


E (pJi + E"”" ' Po-E 

The expression (|5.16|) is extended to the other processes by using Table 1. We have 


reported the results in Figure 8. 

5.3 Photon Emission and Absorption 

In order to eliminate the infrared divergences due to the radiative diagrams of Figures 
5(a) and 5(b), it is necessary to include the rates of processes in which a photon is either 
absorbed or emitted. In fact, they exactly provide the divergent contributions which 
cancel the corresponding ones due to vertex and wave function renormalizations. The 
processes give a hnite contribution as well, which takes the following form|^ 

GUCl + 3Cl)a Ip''' 


+ n ^ e + p) = 


271^ 


TT Jo 


poo poo T-\' 

/ d\p\ / d|k| ^ 5(|k|) 

Jo Jo Pq 
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/ 2 / f \ \ / 

J’o _log<P" + IP 




+ 


-wlog 


PT 

p'o + Ip' 


Q+ + Q'-2Q2] (5.17) 

Ql + Q-1 1, 


2|p'r'° vpo - ip'I. 

where Q\ = (go±|k|)^E,(go±|k|) 0(go±|k|), = qlEy{qo) 0(go) and S = 1. The 

expression (|5.17|) is extended to the other processes by using Table 1 and with S' = 1 for 
Ve + n ^ e~ + p and n e~ Ve + p, and S = 0 for + n ^ Vf. p. The ratio 
Auj^/{a/Ti)u}B is plotted for all processes ([1.2|) in Figure 9. The results agree with the 


ones reported in Figure 8 of Ref. ||I^ apart from the channel e“ -|- p — >■ W n, which 
is there underestimated by a factor 10. Notice that Auj^ represents the leading thermal 
correction to neutron decay and the inverse reaction, in the temperature range relevant 
for BBN. This is simply understood since the inclusion of the process ''y + n ^ e~ + Vf. +p 
greatly increases the neutron decay rates, otherwise strongly suppressed by phase space. 


®The result for Awy quoted in [Q as Eq. (11) has a missing factor 1/E. 
®Our result slightly differs from Eq. (13) in Ref. fl^. 
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6 Analysis of the results and ^He abundance 

In the previous sections we have reported a detailed analysis of three kind of corrections 
to Born rates for n ^ p processes: zero temperature radiative corrections, hnite nucleon 
mass and thermal radiative effects. With reference to Figures (3), (4) and (6)-(9) the 
results for each reaction channel can be summarized as follows. 

(a) z/e + n —e“ + p 

For the crucial BBN temperature range, 0.1 MeV < T < 3.5 MeV, the two main 
corrections to the Born rate come from zero temperature radiative and kinetic terms. 
The contribution Aur is weakly depending on T and represents the dominant term, 
though for large temperature the kinetic correction Auk starts contributing signif¬ 
icantly. The two combined contributions correct the Born rate for a factor 6 -i- 9%, 
whilst thermal radiative ones are of the order of 1%. 

(b) e~ + p ^ Ue + n 

For this channel radiative corrections are dominant at low temperature, while the 
kinetic ones give a quite relevant effect in whole BBN temperature range and correct 
the Born rate for a factor varying from 1% at low T up to 3% for T = 3 -i- 4 MeV. 
The radiative contribution is quite rapidly decreasing with temperature, reaching 
large negative values. This leads to a partial cancellation between Aojr and Aojk- 
Thermal corrections are dominated by bremsstrahlung contribution Auj^ and can 
be as large as 2% of uir for large temperature. 

(c) e'^ + n ^ Ve + p 

For this process the radiative corrections have a behaviour quite similar to chan¬ 
nel (6), though they are even more rapidly decreasing with temperature. This 
again leads to a partial cancellation between Aujr and the positive monotonically 
increasing Aujk- Thermal corrections are again mainly provided by photon emis¬ 
sion/absorption and monotonically increase with temperature up to a factor 2% of 
the corresponding Born rate. 
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(d) z/e + p —> e+ + n 

The radiative and kinetic corrections snm np to a factor of abont 5% of the Born 
rate. Actnally the opposite behavionr of Aur and Auk conspires to give an al¬ 
most constant total correction in the whole interesting temperature range. Thermal 
effects are quite small, contributing for less than 0.5%. 

(e) n ^ e~ +Ve + p 

For the neutron decay the radiative corrections are practically constant and give 
the leading effect for small temperature, T < me- For larger T the thermal photon 
emission/absorption processes rapidly become dominating over all other correction 
terms with a relative ratio to the Born rate as large as 10"^. However, this large 
correction is weakly contributing to total n ^ p rate, and so to abundance 
prediction, since in the temperature range T > me the scattering processes {a)-{d) 
are largely dominant over decay and inverse decay. 

(f) e~ + Ve + p ^ n 

Same considerations of the direct process (e) hold for inverse neutron decay. As 
for the direct channel, the kinetic contribution is negative down to temperatures of 
the order T ~ 0.2 MeV. The contribution of lAcUi^-l is however negligible for both 
processes, smaller than 1%. 

All results are summarized in Figure 10, where we have shown the total relative cor¬ 
rections in percent. For the processes (a) and {d) the correction is almost constant over 
the entire considered range for T, and of the order of 6 -i- 10% and 5 -i- 6%, respectively. 
The positive kinetic contribution soften the deep decreasing of the radiative terms for 
channels {b) and (c). Finally the large effect of thermal bremsstrahlung for neutron decay 
and inverse process (e) and (/) is particularly evident. 

We have shown the total rates for n ^ pin Figure 11, while the total relative correction 
Au/ub = in percent, are plotted in Figure 12. For future applications in the 

BBN codes we have performed a fit of the numerical results for uj{n ^ p) and uj{p ^ n). 
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The fitting expressions are the following 

1 ® ^ T V 

—’ ( 6 - 1 ) 
^ /=0 \ifle' 

^(P^n) = ^exp(-i^) g6,(T) . (6.2) 

where for n ^ p 

ao = l ; ai = 1.0988 ; oa = -9.8297 ; as = 17.379 ; 04 = 32.197 ; 

05 = 27.372 ; og = -0.20975 ; 07 = 0.67033-10-2 ; og = -0.91503-10-^ , (6.3) 

while for p —> n 

61 = 22.138 ; 62 = -77.737 ; 63 = 130.76 ; h = -22.010 ; 65 = 50.262 ; 
be = -4.8032 ; 67 = 0.53087 ; bg = -0.03403 ; 69 = 0.11508-10-2 ; 

610 = -0.15850-10-^ ; q = 2.9123 . (6.4) 


The ht has been obtained requiring that the htting functions differ by less than 0.1% from 
the numerical values. 

The whole correction Au for n ^ p turn out to be a positive decreasing function over 
the whole temperature range relevant for BBN. The main contribution at low temperature 
for both total rates comes from the radiative corrections, while for T > 2 3 MeV kinetic 

contribution starts dominating. This is particularly evident by looking at Figures 13 and 
14, where Aojr/ojb and Aujk/ojb are separately plotted. While for radiative corrections 
the effect on n —> p total rate is larger than on the p —> n one, Aujk shows an opposite 
behaviour. The competition of these two corrections is then responsible for the presence 
of the inversion point in Figure 12 at T ~ 0.15 MeV. Finally, the pure thermal radiative 
corrections 


Aujj' = Aum T Aujw + Auy + Au^ 


( 6 . 6 ) 


are plotted in Figure 15. Their order of magnitude is sensibly smaller, but they never¬ 
theless may contribute for a factor 0.2 -h 0.4% at the freeze-out temperature T ~ 1 MeV. 
Although some differences in the expressions for several contributions, our result for Auj^ 
essentially agree with the results of Ref.s [0,114 . 
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The corrections on the n ^ p rates, as well as the plasma electron mass correction 
to neutrino temperature, considered in Appendix A, may produce a sensible correction 
to the mass fraction 1^, which is strongly dependent on the neutron fraction = 
nn/{nn + Up) at the nucleon freeze-out. An accurate theoretical prediction for Helium 
abundance, as well as for the other light nuclei produced during BBN, can be obtained 
by using the standard BBN code 0 suitably modihed to take into account all the above 


corrections Here we report the results for the expected variation of the surviving 
neutron fraction induced by the whole effects Auj{n ^ p). This allows for a simple 
estimate of the corresponding variation of Yp. 

The neutron fraction is determined by the differential equation 


^ ^ \.^(P Xr,) - uj(n p)Xr,] 


( 6 . 6 ) 


Writing A„ = X^ + 6Xn, where 6Xn is the correction induced by Au, we have at hrst 
order 


-^SXn = -^[{ub{p ^ n)+UB{n ^ p))5Xn-Au{p ^ n) 
dl dl 

+ {Au}{p ^ n) + Aoj{n ^ p)) Xl^ 


(6.7) 


Notice that the zero order abundance has been dehned as the one obtained by the 
Born amplitudes rescaled by the constant factor 961/886.7, which provides at tree level 
the correct prediction for neutron lifetime (see our discussion in section 3). The Born 
rates ujb{p n) in ( |6.7D are therefore rescaled by the same factor. Equations (|6.7|) have 
been numerically solved using our htting function for uj{n •«-> p) and a similar one for 
LOsin ^ p), which we do not report for brevity. We found for the asymptotic abundance 
6Xn — 0.0024, with a relative change, in percent, dX^/X^ = 1.6%. It is also interesting 
to evaluate the effect due to the thermal radiative contribution. Using again (|6.7|) , with 
now the AujT{n ^ p) in the r.h.s., we get = —0.0002, 5X'^/X^ = —0.15%. 

These results allows for a simple estimation of corrections to ^He mass fraction P5 

tn. 


SYp = 2 5Xn exp 


exp 
I n 


(6.8) 


26 

















where tns — 180 sec corresponds to the onset of nucleosynthesis and the exponential factor 
accounts for the depletion of relic neutrons at freeze-out due to /3-decay. Using the results 
for 6X„ we hnd 


SYp 

~ —0.0004 , 

SYp 

yp 

= -0.15% , 

thermal 

SYp 

~ 0.004 , 

SYp 

yp 

= 1.6% , 

total 


(6.9) 

The thermal radiative corrections decrease the value of Yp, in agreement with the results 
of [I^l and |M. However, the total correction is largely dominated by zero temperature 
radiative and hnite mass corrections, which give a positive contribution to Yp. The total 
correction, evaluated with the above method, agrees with the recent result obtained in 
3^ . They also include the effect of hnite temperature QED corrections to the equation 


of state of the electromagnetic plasma and of residual coupling of neutrinos to the 
thermal bath during the annihilation phase |]^, |^. However, the two combined 

corrections give a very small positive contribution, of the order 6Yp ~ (2 3)T0“^. 

7 Conclusions and Outlook 

In this paper we have performed a complete analysis of the corrections to the tree level 
rates ( [1.2|) . These reactions are relevant to hx the neutron density at the freeze-out, 
which is the hrst step of the primordial synthesis of light nuclei. In particular we have 
considered order a radiative corrections, 1/M/v effects obtained relaxing the assumption 
of inhnite massive nucleons, and thermal radiative effects on the rates as well as on the 
behaviour of neutrino temperature after decoupling. The total correction Au to the rate 
per nucleon is a monotonically decreasing function of the temperature in the interval 
0.1 MeV < T < 3.5 MeV, for both processes n <r-> p. It varies in the range (2 -h 8)% of 
the Born rates. The main effects come from radiative and kinematical terms Aujr and 
Aujk, considered in section 3 and 4, respectively. We find that Aujr contributes for about 
(—1-=-8)% in the above range for T, while Auk, which originates from the relative motion 
of the initial nucleon with respect to the comoving frame, grows with temperature up to 
a value of about 2.5% of oosin •«-*• p). Thermal radiative corrections Aut, evaluated in 
the real time formalism, rapidly increase with temperature and only contribute for about 


27 

































(0.1 ^0.7)%. 


We have calculated how these corrections affect the theoretical prediction for '^He mass 
fraction Yp. We have solved the differential equation governing the evolution of neutron 
fraction Xn and found for the total correction induced on Yp by Acn the value SYp ~ 0.004, 
which corresponds to a relative effect 5Yp/Yp ~ 1.6%. This result has been obtained using 


the simple method outlined in |^. The thermal radiative effects alone contribute to this 
variation for a negative factor 6Yp = —0.0004, in agreement with the results of and 
T^. Finally, the total correction to ^He abundance is in substantial agreement with a 


recent estimate El , though the result for the thermal radiative small correction is of the 
same magnitude but with an opposite sign. A more careful analysis of the effect of Ao; 
on Helium and all other nuclei abundances, using a modification of the BBN codes 0 , is 
in progress . 
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A Neutrino temperature versus photon temperature 


In this Appendix we evalnate the ratio of nentrino temperatnre Ty over the photon tem- 
peratnre T. 

When weak interactions go ont of eqnilibrium, at T ~ 1 MeV, neutrinos decouple 
from the photon and thermal bath. Their temperature then decreases as 

with a{t) the cosmological scale factor. Shortly after, as T drops below the electron mass, 
electrons and positrons annihilate and transfer their entropy to photons only. To evaluate 
the ratio Ty/T we recall the expressions for the entropy density of a particle specie i, with 
Qi internal degrees of freedom 


9i rn3 n't^ \ — 9i rp3 


OTT^ OTT^ 


^ (y^ + 2yxiy ^ (Ay^ + 8yxi + ?>xfj (exp(a:i + ?/)±l) ^ dy , 

(A.l) 


where +, — refers to fermions and bosons, respectively, and Xi = rrii/Ti. Using the fact 
that the entropy per comoving volume is constant, and in particular that neutrino entropy 
remains unchanged after decoupling one gets 


Ty / 27r"^ + 15 X(xe) 
Y ^ i^27r4 + 15 X(a:f) 


1/3 


with 


Xe = 


1 ( 

IJle H- 


X? = 




X 


me + 


D 


an ' 

me 


(A.2) 


(A.3) 


The electron neutrino decoupling temperature, evaluated in |^, is To — 2.3 MeV. We 
have verihed, that changing To in the range 2 3 MeV the ratio Ty/T only changes for 

a factor less than 0.2%. We have also included the effect of thermal electron mass on 
neutrino temperature via an effective simple estimate of the form aX^/me . 


B Finite nucleon mass corrections to weak ampli¬ 
tudes 


In this Appendix we compute the spin-summed squared amplitudes for processes (|L2|) for 
finite nucleon masses. 

Let us consider the most general four-fermion effective interaction, responsible for the 
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process Pi + P 2 ^ Ps + P4, 


M = ^u{p4)Oi,u{p2)u{p3)Y{^-'JbMpi) , (B.l) 

where the effective weak coupling is given by 

= 7/.(Cv - C'aTs) + i^(Ti,v{p4-P2T , (B.2) 

Mn 

with = (i/ 2 ) [ 7 ^, 7 ,^]. Adopting the same notation of section 2 we obtain 

y/ |M|^ ~ + 2m\m\ + mlml + 2 7712^4 + 2m 27774 

spins 

+ m\m\ + 2 ml ml + mlml — 2 mis — 2mls — 2mls — 2 mis 
+ 2s‘^ — mlt — mlt — mlt — 2m2m4t — mlt + 2st + t^) 

+ Cv'^ {ml ml + 2 ml ml + ml ml — 2 ml m 2 m^ — 2 m 2 ml m^ 

+ mlml + 2 ml ml + mlml — 2 mis — 2mls — 2mls — 2 mis 
+ 2s^ — m^t — m^t — mgt + 2m2m4t — m^t + 2st + 

2 / 27 22 22 22 2 2 

— 2 Ca - rT ~ ^2 ~ ■'^2 ^3 + ^2 ^4 ~ ^2 ^3 ^4 ~ ^2 7^4 

Mm 

+ m2 mg m4 — m^ m4 + mg m4 — m^ m 2 1 — m^ t 

— m 2 ml t — ml m4 1 — ml m^ t — ml m^ t — m 2 ml t 

— ml t + 2 m 2 st + 2 m 4 st + m 2 + m^ t^) 

+ 2CvCa{ — ml ml + ml ml + ml ml — ml ml + ml t 
+ mlt + mgt + mlt — 2 st — t^) 

, n f‘2 f 2 z , 2 2 4 4 22 

+ 2 Cv -TY~ ^2 + ^2 TUg — m2 mg — m^ m^ — m^ m2 m 4 

Mm 

+ m^ mg m4 — m 2 mg m^ + mg ml — ml m 2 s + m 2 ml s 
+ ml m 4 s — ml m^ s — ml m 2 t — ml t + ml m^ t 

— mlm^t + m 2 mlt — mlt + m 2 t^ + m^t^)). (B.3) 

In the above formula we have used the following dehnition of the two independent Man¬ 
delstam invariants 

s = (pi + P2f = {P3 + Pif , 

t = {pi - pzf = {p2 - Pif . 
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(B.4) 

(B.5) 



In the formula ( |B.3| ) we have retained only the leading terms, up to hrst order in the 
coupling / 2 . 

For the process + n —»• e“ + p, the corresponding X^spins \M\^ is obtained from (|B.3|) 
by replacing mi, m 2 , m 3 and m 4 with 0, M„, me and Mp, respectively. Furthermore, from 
the dehnition of the Mandelstam invariants and using the four-momentum conservation 
we have in the neutron rest frame 


+ 2M„go 5 

(B.6) 

Mp — + 2M„(pg — go) • 

(B.7) 


By using ( |B.3|) and changing accordingly the above correspondences one easily obtains 
the expression for Z^spms for all other processes in (|1.2|) . 


C Kinematical and thermal mass shift corrections to 
the rates 


In this Appendix we generalize to all other processes in ( p..2|) the results reported in (|4.17|) 
and ( b.l4|) for kinematical and thermal mass corrections to Born rates for i>e ^ e~ +p. 

The kinematical corrections of section 4.2 can be obtained from ( |4.17|) by using the 
substitution rule of Table 1 and replacing the factor in square brackets as follows. First 
of all, note that for the reaction z7e + p —»■ -|- n this factor is the same as in (^4.17|) . 
Differently, for e~ + p ^ Ue + n and e'^ + n ^ Vf. + p, one has 




-3-^F,(go) [l-2F,(go)] + 


IP 


rp2 


[1 - nW)] |1 - iFM)] + 3^F,(qo) 

-L IJ 




^ O I / I 

2|p I go 


C2+3C1 


STT^p’q 


2 FM [1 - FM)] 


(C.l) 


while for n 


e Fi>e +P 


4 


-3 - [1 - 2f;(,„)i - n(p;) |i - 2 f ,( p '„)\ + 3 ^fm 


go 


T [c^ + 3Cl) [sTT^p'o, 


^e(Po)i^i.(go) 


(C.2) 


Finally for e + Ug + p 


n 


go 


-3 + ^ |1 - [1 - 2FM] + 


Ip 


'|2 


[1 - Fe(pi)l |1 - 2FM)] - 3^ [1 - F.(9o)] 


2^2 


T 

u 
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^ o I f\^ 

2|p I go 


h - F - f 

T ^ 1^ Q + 3^2 y "iTT^v'^ 


2 [1-F,(p')] [l-F,(go)] 


(C.3) 


A similar replacement is required to generalize the result of (^.141). In particular, in 


addition to the substitutions of Table 1, the factor in (|5.14|) 


, -^e(Po) _ 1 - Fv[(Iq) ' 
,go T Ty 


(C.4) 


should be substituted in the following way. As for kinematical corrections, this factor 
applies to the reaction Ve +p —*• e^ + n as well. For e“+p —> z/g +n and e+ + n —>• Pg +P, 
one has instead 


while for n —> e + Ug + p 


^2 l-F,{p',) F^iqoY 

,go T n 


^^^e{Po) F^iqoY 


(C.5) 


go 


T 


T 


(C.6) 


and e + z/g + p —> n 


^ _ 1 - FejPo) 1 - F^iqoY 

go T T, 


(C.7) 


References 


[1] P.J.E. Peebles, Physical Cosmology, (Princeton University Press, 1971, Princeton); 

S. Weinberg, Gravitation and Cosmology, (Wiley, 1972, New York); 

D.N. Schramm and R.V. Wagoner, Ann. Rev. Nucl. Sci. 27 (1977) 37. 

[2] E.W. Kolb and M.S. Turner, The Early Universe, (Addison-Wesley Publishing Com¬ 
pany, 1990, New York). 

[3] K.A. Olive, Proc. of the Workshop TAUP97, Gran Sasso Laboratory, Italy, September 
7-11, 1997, |astro-ph/971T[60| ; 

G. Steigman, |astro-ph/9803055 . 

[4] K.A. Olive and D. Thomas, Astropart. Phys. 7 (1997) 27; 

Y.I. Izotov, T.X. Thuan, and V.A. Lipovetsky, Ap. J. Suppl. 108 (1997) 1; 

R.F. Garswell, M. Rauch, R.J. Weymann, A.J. Gooke and J.K. Webb, MNRAS 268 


32 




















(1994) LI; 

A. Songaila, L.L. Cowie, C. Hogan and M. Rngers, Nature 368 (1994) 599; 

S. Bnrles and D. Tytler, Proc. of the ISSI Workshop Primordial Nuclei and their Galac¬ 
tic Evolution, |astro-ph/9712265 . 

[5] S. Sarkar, Kept. Prog. Phys. 59 (1996) 1493. 

[6] L.M. Kranss and P. Romanelli, Ap. J. 358 (1990) 47; 

P.J. Kernan and L.M. Kranss, Phys. Rev. Lett. 72 (1994) 3309; 

G. Fiorentini, E. Lisi, S. Sarkar and F.L. Villante, |astro-ph/9803177 . 

[7] R.V. Wagoner, Ap. J. Suppl. 18 (1969) 247; Ap. J. 179 (1973) 343; 

L. Kawano, preprint FERMILAB-Pnb-88/34-A; preprint FERMILAB-Pnb-92/04-A. 

[8] A. Sirlin, Phys. Rev. 164 (1967) 1767. 

[9] W.J. Marciano and A. Sirlin, Phys. Rev. Lett. 56 (1986) 22, and references therein. 

[10] W.J. Marciano and A. Sirlin, Phys. Rev. Lett. 46 (1981) 163. 

[11] D.M. Wilkinson, Nucl. Phys. A377 (1982) 474. 

[12] C. Caso et ah, Eur. Phys. Jour. C3 (1998) 1. 

[13] D.A. Dicns, E.W. Kolb, A.M. Gleeson, E.G.G. Sndarshan, V.L. Teplitz and M.S. 
Tnrner, Phys. Rev. D26 (1982) 2694. 

[14] J.L. Gambler, J.R. Primack and M. Sher, Nucl. Phys. B209 (1982) 372. 

[15] L. Dolan and R. Jackiw, Phys. Rev. D9 (1974) 3320. 

[16] J.F. Donoghne, B.R. Holstein and R.W. Robinett, Ann. Phys. (N.Y.) 164 (1985) 23. 

[17] J.F. Donoghne and B.R. Holstein, Phys. Rev. D28 (1983) 340; Phys. Rev. D29 (1984) 
3004. 

[18] A.E. Johansson, G. Peresntti and B.S. Skagerstam, Nucl. Phys. B278 (1986) 324. 

[19] W. Keil, Phys. Rev. D40 (1989)1176. 


33 








[20] R. Baier, E. Pilon, B. Pire and D. Schiff, Nucl. Phys. B336 (1990) 157. 

[21] W. Keil and R. Kobes, Physica A158 (1989) 47. 

[22] M. LeBellac and D. Poizat, Z. Phys. C47 (1990) 125. 

[23] T. Altherr and P.Aurenche, Phys. Rev. D40 (1989) 4171. 

[24] R.L. Kobes and G.W. Semeneff, Nucl. Phys. B260 (1985) 714; ibidem B272 (1986) 
329. 


[25] R.F. Sawyer, Phys. Rev. D53 (1996) 4232. 

[26] LA. Chapman, Phys. Rev. D55 (1997) 6287. 

[27] S. Esposito, G. Mangano, G. Miele and O. Pisanti, Phys. Rev. D58 (1998) 105023. 

[28] M.A. Herrera and S. Hacyan. Ap. J. 336 (1989) 539; 

N.C. Rana and B. Mitra, Phys. Rev. D44 (1991) 393; 

S. Dodelson and M.S. Turner, Phys. Rev. D46 (1992) 3372; 

A.D. Dolgov and M. Fukugita, Phys. Rev. D46 (1992) 5378; 

N.Y. Gnedin and O.Y. Gnedin, |astro-ph/9712199| ; 

S. Hannestad and J. Madsen, Phys. Rev. D52 (1995) 1764; 

A. D. Dolgov, S.H. Hansen and D.V. Semikoz, Nucl. Phys. B503 (1997) 426; 

B. Fields, S. Dodelson and M.S. Turner, Phys. Rev. D47 (1993) 4309. 

[29] D. Seckel, preprint BA-93-16, |hep-ph/9305311| ; 

R. E. Lopez, M. S. Turner and G. Gyuk, Phys. Rev. D56 (1997) 3191. 


[30] V.P. Gudkov and K. Kubodera, |nucl-th/9706074 . 


[31] F. Halzen and A.D. Martin, Quarks and Leptons, (Wiley, 1984, New York). 

[32] S. Esposito, G. Mangano, G. Miele and O. Pisanti, in preparation. 


J.Bernstein, L.S. Brown and G. Feinberg, Rev. Mod. Phys. 61 (1989) 25. 


[34] R.F. Lopez and M.S. Turner, preprint FERMILAB-Pub-98/232-A (1998), |astro- 
I ph/9807279[ 


34 












[35] A.J. Heckler, Phys. Rev. D49 (1994) 611. 

[36] K. Enqvist, K. Kainulainen and V. Semikoz, Nucl. Phys. B374 (1992) 392. 


35 




36 



sec 



Figure 2: The Born rates ujb for all n ^ p processes 
correspond to the reactions as labelled in (p..2D . 


(see sect. 2). Hereafter the curves 
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gure 13: ihe zero temperature radiative correc 
ashed line) processes. Finite mass contributioi 
id weak magnetism are also included. The resu! 










